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Abstract 

Both amorphous and crystalline materials frequently exhibit low temperature specific heats in 
excess of what is predicted using the Debye model. The signature of this excess specific heat is 
a peak observed in C /T^ versus T. To understand the curious absence of long-range ordering of 
local distortions in the crystal structure of pyrochlore Bi2Ti207, we have measured the specific 
heat of crystalline Bi2Ti207 and related compounds. We find that the peak in C/T^ versus T 
in Bi2Ti207 falls at a substantially lower temperature than other similar compounds, consistent 
with the presence of disorder. This thermodynamic evidence for disorder in crystalline Bi2Ti207 
is consistent with quenched configurational disorder among Bi lone pairs produced by geometrical 
frustration, which could represent a possible realization of "charge ice" . 
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It is well known that at relatively low temperatures, typically in the range of 2 K to 
30 K, amorphous systems exhibit a larger specific heat than the simple prediction of the 
Debye model. This excess specific heat manifests itself as a peak in C /T^ versus T, and 
is generally attributed to local low energy vibrational modes not accounted for in the De- 
bye model. These low energy modes are observed in Raman spectrap] and from inelastic 
neutron scattering, [2] where they give rise, in glasses, to the so-called "boson peak".[3l H] 
The mechanisms responsible for these low energy modes have been discussed in terms of 
localized vibrations ("floppy modes"), domain wall motions of the glassy mosaic structure, 
and transverse phonon modes. |S1 El |7| 

A number of crystalline materials also exhibit a low-temperature peak in C/T^. This peak 
can be attributed to van Hove singularities where the vibrational density of states (VDOS) 
crosses the Debye density of states, leading to a flattening of phonon dispersion curve. [8] 
This is nearly equivalent to stating that the local modes are responsible; in materials such 
as Si02, the nature of the vibrations responsible for the peak in crystalline and amorphous 
samples are similar. [3] As low temperature specific heat measurements can provide evidence 
for disorder in insulating crystalline compounds, we have used thermodynamic measurements 
to investigate displacive disorder in Bi2Ti207 and related systems. 

In this contribution, we examine the excess low temperature specific heat in a series 
of structurally- and compositionally-related crystalline compounds of interest as polar and 
dielectric materials. The pyrochlore compound Bi2Ti2060' is compositionally and electron- 
ically related to the perovskite ferroelectric PbTiOs, in that the A site of the structure is 
occupied by a main group ion with the 6s^, configuration, and the B site is occupied by 
dP Ti^"*" (structures displayed in FIG.[T|. Both the A and the B site ions in these compounds 
are therefore susceptible to off-centering. This so-called stereochemical activity, and this is 
what drives the phase transition to a coherent, polar, tetragonal ground state in PbTiOa. 
The related pyrochlore compound Bi2Ti2060' (more simply referred to as Bi2Ti207) shows 
no such phase transition and even at 2K, the crystal structure is cubic. [9] However, in a 
manner that is common across many pyrochlore crystal structures with Bi^"*" on the A site, 
Rietveld refinement of neutron diffraction data indicates that Bi is locally off-centered. This 
evidence for local distortion, in conjuction with the absence of ferroelectricity, suggests that 
that the Bi ions are displaced randomly into one of several distinct, but symmetric, sites 
about the equilibrium position as indicated in FIG.[T](b).[ni [IHl El ESI [ISl IS] 
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FIG. 1: (a) Structure of ideal Bi2Ti2060'. The lone pair around Bi'^+ atoms are visualized using 
the valence electron localization function (ELF) with isosurfaces of value 0.625 following methods 
presented in reference [151 (b) Split atom positions of Bi and O' (50% isotropic thermal ellipsoids 
displayed) as obtained from Rietveld refinement of powder neutron diffraction data indicative of 
local displacive disorder, following reference [TOl (c) Crystal structure of tetragonal PbTiOs (two 
unit cells displayed), using the ELF to locate lone pairs around Pb^"*". Note the difference in the 
nature of the lone pairs in (a) and (c). Black spheres are Bi or Pb, cyan spheres are Ti, and orange 
spheres are O. 

The goal of this work is to explore whether there is any thermodynamic evidence for dis- 
placive disorder in Bi2Ti207 driven by geometric frustration on the pyrochlore lattice, [16] 
Frustration of ferromagnetic Ising spins in compounds with the pyrochlore crystal structure 
are well known. [T7[ [TS] with analogies having being drawn with the problem of hydrogen 
positions in the crystal structure of cubic ice Jc, with its associated residual entropy. [I9] A 
similar analogy between ice and putative polar pyrochlores was suggested in reference [151 
that Bi2Ti207 may be a manifestation of "charge ice". We emphasize that this disorder 
is expected to be driven by purely geometrical considerations and develop on well-ordered 
crystalline lattices, so this "charge ice" structure should be distinguished from the "charge 
glass" state observed in structurally disordered La2Cui„^.Li2,.04 and Lao-^Sr^NiO/i . [20] 

In order investigate the role of disorder on the low temperature heat capacity in crys- 
talline systems we have measured the specific heat of Bi2Ti207, two isostructural compounds 
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Y2Ti207 and Bi2NbIn07, as well as the ferroelectrics PbTiOs and Bi4Ti30i2. Coherent lone 
pair displacements in the latter two systems produce a ferroelectric ground state, which sug- 
gests that residual displacive disorder should be minimal. The compound Y2Ti207 unlike 
the four others, has no lone pairs on the A site and can be expected to display the least dis- 
order. All the samples are white to light yellow and highly electrically insulating in powder 
form. 

The polycrystalline samples studied here, with the exception of Bi2Ti207, were prepared 
by ceramic routes and phase purity was verified using powder X-ray diffraction. Perovskite 
PbTiOs was prepared from PbO and Ti02, with a 5% stoichiometric excess of PbO to balance 
losses due to volatilization. Pyrochlore Y2Ti207 was prepared from Y2O3 and Ti02. The 
Aurivilius phase Bi4Ti30i2 was obtained from appropriate stoichiometric starting ratios of 
Bi203 and Ti02. Bi2lnNb07 was prepared by reacting appropriate stoichiometric amounts 
of Bi203, Nb205, and In203, following the method reported by Zhou et a/. [13] The Bi2Ti207 
sample used in this study was prepared by low temperature routes and previously subject 
to neutron structure determination. [H] All specific heat measurements were carried out on 
Quantum Design Physical Properties Measurement Systems (PPMS). We used different 
techniques to extract the low temperature specific heat of these compounds. The Bi2Ti207, 
PbTi03, and Bi4Ti30i2 powders were pressed into dense pellets, mounted to a calorimeter 
with thermal grease. The Bi2NbIn07 powder was cold sintered with Ag powder, with the 
specific heat contribution from the silver measured separately and subtracted. The specific 
heat of Y2Ti207 was measured both by mixing the loose powder with a small amount of 
thermal grease then compacted between aluminum plates with weighing paper and also by 
pressing it in a 1:1 ratio with silver. Both techniques yielded quantitively similar results, 
so we present only the data from Y2Ti207 mixed with thermal grease. The Bi2Ti207 and 
Y2Ti207 samples were separately measured down to 0.5 K using a PPMS system equipped 
with a ^He refrigerator, with the relaxation measured over two time- const ants. 

The specific heat capacities for all five samples are shown in FIG.[2|a). In order to more 
clearly emphasize the low energy features in the system, giving rise to specific heat in excess 
of the Debye contribution, we have plotted the specific heat as C /T^ versus T. All samples 
show a maximum when plotted in this manner, with the amplitude of the peak increasing 
monotonically with decreasing peak temperature. This peak indicates that the specific heat 
of these crystalline samples exceeds the Debye C(T) at low temperatures, and are consistent 
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FIG. 2: (color on-line) (a) Temperature dependence of the low temperature specific heat (plotted 
as C /T^ versus T) for the five samples studied in this contribution. The values are per mole of 
the A-site atoms in the different formulae. For all samples, the values of C /T^ in excess of about 
0.5 mJ mole~^ are contributions not accounted for by the Debye model, (b) Scaled plots: 
C/T^/Pc as a function of the scaled temperature T/Tm- 

with previous measurements on crystalline PbTiOa, which has been shown to exhibit a peak 
at Tm = 15K.|2l] 

As Y2Ti207 and Bi2Ti207, have the same crystal structure, one would expect that both 
compounds should exhibit similar vibrational DOS. Contrary to this expectation, we find 
that Bi2Ti207 has a much larger low energy density of states, which is partly due to the larger 
mass of Bi compared to Y. This larger mass cannot, however, fully explain the appearance of 
a low-energy peak; we attribute this excess heat capacity to the presence of disorder among 
Bi lone pairs. The lower temperature for the Bi2Ti207 peak relative to the Y2Ti207 peak is 
consistent with increasing disorder in the Bi pyrochlore arising from lone pair displacements, 
as discussed in more detail in the following. The other pyrochlore compound studied here, 
which is also expeced to displacive disorder on the A site, Bi2NbIn07, shows a relatively 
large, low-temperature peak in C/T^, although the amplitude of the peak, Pq is smaller, 
and the temperature Tm at which the peak maximum is found is higher that Bi2Ti207. 
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In order to compare the excess specific heat in the two pyrochlore compounds with other 
crystalhne lone-pair active systems, we have also measured the low temperature specific 
heat of PbTiOs and Bi4Ti30i2. It is seen in FIG.[2]^a) that these compounds also display a 
distinct peak in C/T^, although the peaks fall at higher energy, and are smaller in magni- 
tude than the peaks displayed by the two pyrochlore compounds with displacive disorder. 
The distinct difference in the excess specific heat of the lone-pair disordered compounds, 
Bi2Ti207 and Bi2NbIn07 from the lone-pair ordered compounds, PbTiOs and Bi4Ti30i2 
is seen in the scaled C /T^ vs. T plots displayed FIG.[2]^b). The scaling was performed with 
respect to the individual peaks, Pc at which C /T^ is maximum, and the temperature Tm 
where the maximum is found. While materials having both ordered and disordered distor- 
tions exhibit a specific heat peak, the lone-pair disordered compounds display a distinctly 
greater width in the distribution of the excess specific heat. This increase in the scaled full 
width half maximum (FWHM) of the Bi2Ti207 and Bi2NbIn07 curves indicates that the 
modes contributing to the excess entropy in these materials have a wider distributions of 
energies than in the ordered compounds. As a reference, the curve computed for a single 
mode Einstein oscillator is plotted as a solid line in FIG.|2]^b). The relative width of this 
excess heat capacity peak, plotted on a semilog scale and measured at 70% of the maximum 
value, increases from 0.36 for PbTiOa to 0.71 for Bi2NbIn07 and 1.05 for Bi2Ti207. This 
increase is similar to the broadening of the C/T^ maximum in amorphous Si02 relative to 
crystalline 8102^3' ^^"^ t)e associated with increasing disorder of the lone pair electrons 
in Bi2NbIn07 and Bi2Ti207. 

To quantitatively investigate these peaks, we fit the specific heat in FIG.|3]to the sum 
of a Debye contribution (C/j), with two additional low energy Einstein modes {Ce)- This 
is equivalent to supposing that the VDOS can be modeled by a delta function as it crosses 
the Debye DOS, rather than the finite jump expected for van Hove singularities in 3D 
crystals. We included two separate Einstein modes as FIG.[2|^b) indicates that the excess 
heat capacity cannot be fit by a single mode. While simplistic, this model allows us to relate 
the thermodynamic anomaly to microscopic properties in the samples. For all crystalline 
samples included in this study (Bi2Ti207, Y2Ti207, PbTiOs, Bi2NbIn07, and Bi4Ti30i2) 
we found that a suitable fit could be obtained using a Debye temperature oi Oj:, = 280 K. 
This emphasizes that the underlying lattice heat capacities of these different samples are 
similar, so that the significant differences in the low temperature heat capacities should be 
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FIG. 3: (color on-line) Fits to the specific heat of Bi2Ti207 using a combinations of Debye and 
Einstein modes, as described in the text. Data for Y2Ti207 are seen to be almost completely fit 
by the single Debye mode using the same value of 0/5 used for Bi2Ti207. 

attributed to differences in the lone pair electron behaviour. 

The low temperature C{T) can be modeled as a sum of Debye and Einstein terms with 
each Einstein oscillator giving a contribution: 



to the specific heat. In this expression, p is the spectral weight, R the gas constant, and houo 
the mode energy. For Bi2Ti207 the first excitation has an energy of Ei = 2.3 meV and an 
oscillator strength of pi = 0.32, while the second excitation falls at E2 = 5.5 meV and has 
an oscillator strength of p2 = 2. Similar fits were obtained for the other samples. We find 
that Bi2Ti207 has the lowest energy oscillator. The values of the oscillator strengths and 
energies for the different samples are given in Tablejlj We find that Hooq for the lowest energy 
mode in these fits depends linearly on the temperature of the peak. Motivated by this result, 
which, within the constraints of our simple model, confirms that the peak temperature and 
vibrational mode frequencies are linearly related, we will discuss our data solely in terms 
of the peak temperature. This assumption allows us to avoid specific details of the fitting 
parameters, allowing a more general and robust analysis of the results. 

To more clearly demonstrate the influence of increasing disorder on the C/T^ peak in 



Ce=pR 



(1) 
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El (meV) 


Pi 


E2 (meV) 


P2 


Bi2Ti207 


2.3 


0.32 


5.5 


2 


Bi4Ti30i2 


3.6 


0.7 


6.7 


5 


BisNblnOy 


2.8 


0.27 


6 


1.6 


PbTiOa 


5.5 


0.4 


8 


1.5 



TABLE I: Energies and number of oscillators used in fitting Einstein modes to the specific heats 
of the diferent compounds. 

heat capacity in crystalline materials, we compare our measurements with some other low 
temperature heat capacity studies on ordered and disordered crystals. We plot the excess 
low temperature heat capacity as C/T^ versus In T for amorphous Si02 and crystalline 
quartz in FIG.|4](a), taken from Ref. [2^, and for Bi4Ti30i2 and Bi2Ti207 in FIG.|4](b). 
The lattice disorder present in amorphous Si02 leads to a reduction in temperature of the 
excess C/T'^ heat capacity peak as compared to structurally well-ordered quartz; this peak 
shift is a signature of increasing disorder in otherwise similar systems. Qualitatively, this 
is similar to the shift in peak temperature between Bi4Ti30i2 and Bi2Ti207, despite the 
fact that both of these materials are crystalline and have very similar Debye temperatures 
and elemental constituents, and hence bare oscillator spectral weights, that are roughly the 
same. These measurements provide thermodynamic justification to support the suggestion 
that crystalline Bi2Ti207 has additional disorder as compared to Bi4Ti30i2. 

Because there is no evidence for any lattice disorder in Bi2Ti207, the presence of a low- 
temperature C/T^ peak provides empirical evidence suggesting that the lone pair electrons 
on Bi are disordered due to geometrical frustration, consistent with previous XRD studies [9]. 
Increasing disorder in other crystalline systems has been observed to shift this heat capacity 
peak to lower temperatures, including oxides with chemical site disorder, such as Ca replac- 
ing Sr in the CaSrFeCoOs brownmillerite[24J. This disorder-induced shift is also observed 
in metallic systems, although the size of the effect may be significantlly smaller than in 
insulators [8j. It has also been noted that geometrical frustration present in systems having 
underconstrained degrees of freedom, such as ZrW208, can also lead to excess low tem- 
perature heat capacity in crystalline systems These results suggest that the downshift 
of the specific heat peak in displacively disordered, but still crystalline, Bi2Ti207 may be 
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FIG. 4: (color on-line) (a) C /T'^ versus In T amorphous Si02 and quartz[23], (c) C/T'^ versus In 
T for Y2Ti207 (squares) and Bi2Ti207 (stars). 

associated with the development of a frozen "charge-ice" state in these systems driven by 
geometrical frustration, rather than the frozen glassy state associated with structural disor- 
der. Furthermore, by investigating the detailed phonon density of states in Bi2Ti207 and 
isostructural Y2Ti207 and correlating these results with heat capacity data, it may be pos- 
sible to gain insight into the microscopic mechanisms giving rise to excess low temperature 
entropy in disordered systems. 

In summary, we find that several insulating crystalline materials with lone pairs, and with 
the pyrochlore, perovskite, and Aurivillius crystal structures, exhibit an excess low temper- 
ature specific heat above the Debye background. The specific heat peaks in pyrochlores 
Bi2Ti207 and Bi2NbIn07 are however, distinctly broader and stronger, and are found to 
appear at lower temperatures than those of isostructural Y2Ti207 and compounds show- 
ing coherent coherent lone pair displacements (PbTiOa and 814X13012). We attribute this 
suppression of the excess specific heat to the presence of additional disorder in Bi2Ti207 
and Bi2NbIn07 arising from incoherent frozen displacements of the Bi lone pair electrons. 
This suggests that Bi2Ti207 may be an example of a structurally well-ordered system with 
charge disorder introduced by geometrical frustration on the pyrochlore lattice. 
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